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Abstract 



We present a status report on our study of long-distance contributions to the decay amplitudes 
A{K^ — > 27r, /) in the framework of the expansion. We argue that a modified prescription 
for the identification of meson momenta in the chiral loop corrections has to be used to gain a self- 
consistent picture which allows an appropriate matching with the short-distance part. Possible 
uncertainties in the analysis of the density-density operators Q% and which dominate the 
CP violation parameter e' je are discussed. As a first result we present the long-distance 
correction to the gluon penguin operator Qg in the chiral limit. 



^Invited talk presented at the 'Workshop on K Physics', Orsay, France, May 30 ~ June 4, 1996. A 
slightly shortened version will be published in the proceedings. 



1 Introduction 



The origin of CP asymmetries is an open issue within the electroweak theory. In the 
standard model CP violation arises from a complex phase introduced by hand in the 
Yukawa couplings. It is of special interest to investigate whether in this approach the ratio 
e'/e, the measure of direct CP violation in K ^ tttt decays, may be close to zero, therefore 
mimicking the superweak model of Wolfenstein [0]. 

At present the experimental evidence is inconclusive, 

p. = / (23 ± 3.6 ± 5.4) ■ 10-^ NA 31 [2] 

[ej \ (7.4 ±5.2 ±2.9) ■ 10-^ E 731 [| ' 

and the new experiments are awaited to provide more accurate data. 

In the standard model the calculation of e'/ e is based on the effective low-energy hamil- 
tonian for AS* = 1 transitions |0, 



= %^«Ec.(/^)Q.(/^), (2) 

v2 i=i 



Ci(/i) = z^ifx) + ryiifx) , r = -^/^n , = V*^^^ , (3) 

where the Wilson coefficient functions Cj(/i) of the local four-fermion operators QiifJ^) are 
obtained by means of the renormalization group equation. They were computed in a 
comprehensive next-to- leading order analysis by two groups ||, |]. Using eq. (0) the ratio 
e'/e reads 

7 = ^mm7^^^*Ei/^((Q«)o - -m2) , 00-' = — ^ 22 . (4) 

{Qi)i are the hadronic matrix elements for the isospin states / which contain the long- 
distance contribution to the amplitudes Aj for the process under consideration, 

{Q,{li))r^{nn,imfi)\K'). (5) 

The dominant CP violating effects result from the gluon and the electroweak penguins, 
(Q6)o and {Qs)2 respectively, with 

= -2 ^(1 + 75)gg(l -75)c? , = -3 J2 egs(l ±75)gg(l -75)^ , (6) 

q=u,d,s q=u,d,s 

where (Q6)2 in addition arises through isospin breaking effects (vrr^r]' mixing). Conse- 
quently, the main task within the investigation of long-distance effects is to estimate the 
degree of cancellation between the two penguin contributions in eq. (^) which gives rise to 
a small value of the ratio e'/e. 

Although recently there has been considerable progress in lattice calculations of the 
hadronic matrix elements 0, it is still meaningful to improve the existing phenomeno- 
logical studies in order to see whether the results are in agreement. There are several 



2 



approximations of non-perturbative QCD available |Ty| . We shall perform our analy- 
sis applying the expansion (A^ being the number of colors) where we will point out that 
a consistent matching between short- and long-distance contributions is possible, provided 
a proper identification of meson momenta in the chiral loop corrections is carried out. 



2 General Framework: l/AT Expansion 

To estimate the hadronic matrix elements we shall use the chiral effective lagrangian for 
K 7T7T decays, which involves an expansion in meson momenta where terms up to 
(9(p^/A^) are included, being the chiral symmetry breaking scale of 0{lGeV) 0, 



X 



(7) 



/ and r are free parameters related to the pion decay constant F.^ and to the quark masses, 
respectively, m = diag{mu, ttLii, rris) being the quark mass matrix. The degrees of freedom 
of the complex matrix U are identified with the pseudoscalar meson nonet, given in a 
non-linear representation. 



U = exp f 7r«Aa 



exp^ 



V2K- 



V2n 

Us 



1^0 



V2K+ 
V2K^ 
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A straightforward bosonization yields the chiral representation of the corresponding quark 
currents and densities, which allows to express the four-fermion operators Qi in terms of 
the meson fields, 



= qa^{l-l^)qj 
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X 



(9) 



The 1/A^ corrections to the matrix elements {Qi)i are calculated by chiral loop diagrams. 
In these diagrams we encounter divergences which are regularized by a finite cutoff Kcut 
as it was introduced by Bardeen, Buras and Gerard (BBG) to estimate the long-distance 
contribution to the A/ = 1/2 rule P]. 

The inclusion of the one-loop corrections leads to a quadratic as well as to a logarithmic 
dependence of the matrix elements on the cutoff, i.e., the results are given in terms of 



A 



cut 



In 



A 2 

^ mi, ' 



(10) 



{j^A,B being pseudoscalar meson masses), where the quadratic part corresponds to the 
chiral limit (m^ = 0). Actually, the loop expansion involves a series in 1//^ ~ 1/^, which 
is in direct accordance to the short- distance expansion in Os/tt ~ 1/A^. 
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Due to the truncation to pseudoscalar mesons, Kcut has to be taken at or even below 
the 0{lGeV). This restriction is a common feature of the phenomenological approaches 
at hand, in which higher resonances are not included. 

To retain the physical amplitudes Aj, which as a matter of principle are scale-inde- 
pendent, the long- and the short-distance contributions are evaluated at the cutoff scale, 
i.e., the long-distance ultraviolet cutoff is identified with the short-distance infrared one. 
This procedure for the matching of the Wilson coefficient functions with the hadronic 
matrix elements we shall analyse in detail in the following section. 



3 Chiral Loop Corrections 

As mentioned above the calculation of chiral loop effects motivated by the 1/N expansion 
was introduced by BBG to explain the A/ =1/2 rule. The authors considered the loop 
corrections to the current- current operators Qi and Q21 whereas the gluon penguin operator 
Qq was included at the tree level [§, implicitly assuming that the 1/A^ corrections to the 
latter are small. 



Following the same lines Buchalla et al. [[TTl performed a detailed analysis of the ratio 



e' /e, in which they stressed that through the quadratic dependence of (Qg) and (Qs) on the 
running mass the evolution of the corresponding Wilson-coefficients is cancelled in the 
large- limit without considering chiral loops. One crucial result of their investigations is 
a small value of e' /e, caused by a large cancellation between the gluon and the electroweak 
penguin contributions. 

Heinrich et al. extended the analysis by including chiral loops, i.e., the 1/A^ corrections 
to the matrix elements of Qg and [jl^ |TB[ . As a net effect they reported an enhancement 
of {Qq)o and a large decrease for {Qs)2, through which the cancellation gets less effective. 



Numerically the authors found a large positive value of e' /e Jll 



e'/e = (9.9 ± 4.1) • 10"^ [m,(lG'e1/) = 175 MeV] . (11) 

The relative size of the 1/A^ corrections contained in eq. (|lTD has to be compared with the 
factors Bq and Bs, denoting deviations from the lowest order results for the corresponding 
matrix elements, as obtained e.g. by lattice calculations. The latter give [0 

5^^=°) = 1.0 ± 0.2 and 5^^=^^ = 1.0 ± 0.2 (12) 

(in accordance with a recent QCD sum rule analysis of {Qq) 0), predicting again a strong 
cancellation between the various penguin contributions. In view of this different results we 
perform a detailed reconsideration of chiral loop effects in the 1/A^ approach. 

One crucial point within our analysis of the K decay amplitudes is the matching 

of short- and long-distance contributions at the cutoff scale. In the standard approach 
referred to above, the effective color singlet gauge boson connecting the two currents and 
densities, respectively, is integrated out from the beginning. The basic diagrams hence to be 
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calculated are presented in fig. 1 where the squares denote weak vertices corresponding to 
the bosonized four-fermion operators Qi and the circles indicate strong interaction vertices. 




Fig. 1: Chiral diagrams for K — > vrvr decays. 

Within this procedure the cutoff Kcut is associated to the meson in the loop, i.e., the 
integration variable is identified with the meson momentum. Consequently, as there is 
no corresponding quantity in the short-distance part, a rigorous matching of long- and 
short- distance contributions is not possible. 

The ambiguity is removed by associating the cutoff to the effective color singlet gauge 



boson as introduced within a study of the K^-Ks mass difference |T5| and recently used by 
Fatelo and Gerard who calculated the long-distance evolution of current-current opera- 
tors in the chiral limit applying the expansion |16[. The corresponding identification 



of the loop integration variable with the momentum flow between the two currents and 
densities, respectively, is feasible in the long- and the short-distance part of the analysis, 
through which a direct momentum matching is justified. 

The crucial feature of this modified approach is illustrated in fig. 2. The momentum of 
the virtual meson is shifted by the external momentum, the former being no longer iden- 
tical to the integration variable g, which affects both, the quadratic and the logarithmic 
behaviour of the 1/N corrections. 




short-distance 




Fig. 2: Matching of short- and long-distance contributions. 

Obviously the modified procedure described above is applicable only for the non-factorizable 
part of the interaction. In fig. 3 we split the first loop diagram of fig. 1 respecting the 
mesonic currents which contribute to the matrix element {'n-^'n-^\Q2\K^) , where Q2 = 

{Jli.)su{JL)ud- 
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+ 




Fig. 3: Representative loop contribution to the matrix element (vr^vr \Q2\K^), the cross circles 
denoting the currents {Jj^)su (left circle) and {Jj^)ud (right) for incoming particles. 

The factorizable loop corrections shown in the third row are completely absorbed through 
the renormalization of the mesonic wave function and the bare decay constant /, a pure 
long-distance feature connected with current conservation and the knowledge of the phys- 
ical decay constant F^. 

Calculating all diagrams of fig. 1, splitted according to a clear identification of the 
virtual meson momenta, the evolution of the matrix element in the chiral limit reads 



7T^n-\Q2iAcut)\K'') = F^ 



m 



K 



2 ^^cut 



(477)2 p 



(13) 



the factor of 3 in front of the cutoff dependent term to be compared with a factor of 2 
obtained in ref. [§]. The result presented in eq. (|T3D is identical to the one derived from a 
study of the operator evolution [|TB|, where the analysis of the matrix element allows us to 
include in addition the logarithmic cutoff dependence arising through non- vanishing quark 
masses (currently being under investigation). 

As illustrated above a rigorous analysis of current-current hadronic operators using the 
modified approach is feasible. The treatment of density-density operators, and 
being dominant for e' /e, requires additional remarks. 

Whereas the factorizable part of the loop diagrams can be removed in the case of 
current-current operators, the corresponding contributions have to be considered explicitly 
for density- density operators; a feature due to the absence of a conservation law for densities 
and, in addition, to the ignorance of a low-energy constant, analogous to the decay constant 

in case of the currents, in which the loop corrections could be absorbed. 

Since a clear identification of the virtual meson momenta from a comparison with the 
short- distance part is not at hand for factorizable diagrams (representing pure long-distance 
phenomena), there is no prescription to fix a possible shift {q ^ q + p)- Consequently, an 
ambiguity in the sub-leading divergent terms occurs. 
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On account of chiral symmetry no quartic dependence on the cutoff is present in the 
matrix elements. Nevertheless, single loop diagrams may produce quartic terms which 
cancel in the sum. Thus a dependence on the (arbitrary) momentum shift may result 
in the (9(A^„J. A concrete analysis however shows that in case of the density-density 
operators each factorizable diagram is at most quadratically divergent, the ambiguity only 
arising in the logarithmic terms. This characteristic is based on the symmetric structure 
of the meson density given in eq. (P), through which the weak vertex solely involves the 
sum of the two (incoming) virtual momenta qi and q2, 



{Dl)^, ^ (d'U),, ^ (^^[Tr'^AJ"),, ^ (gi + g2 + . . -f 



Pext 5 



(14) 



qi + q2 being independent of the integration variable q {n denoting the number of particles, 
Pext the external momentum associated with the density). 

A strict analysis of the non-factorizable contributions can be performed for density-density 
in the same way as for current- current operators. Consequently, from the sum of the 
factorizable and the non-factorizable parts we can clearly determine the quadratic cutoff 
dependence of the matrix elements. 

Splitting the basic diagrams of fig. 1 for the gluon penguin operator Qq = —2{DjiDL)sd 
we may calculate the evolution of the corresponding matrix elements for K ^ titi decays. 
The final result (obtained from about 100 explicit diagrams) reads 



{n+n-\Q,\K') = {it%'\Q,\K') 
where 



-F^(m^ -m^) — 

X 



1 + 3- 



A' 



cut 



f 



1-3 



A 



cut 



(47r)2 /2 
(2m|^ - mlY 



+ Oil/N') , (15) 



(16) 



ml{Aeut) 

Note that (due to the inclusion of the ?7o) the parameter has no quadratic cutoff depen- 
dence. 



It is known that the value of the chiral symmetry breaking parameter A^ increases through 
logarithmic corrections from A^(LO) ^ 1.02GeV to A;^(NLO) ~ 1.20 GeF (see BBG |[), 
thus counteracting the increase caused by the explicit loop contributions given in eq. ([I5|). 
Consequently, we expect the 1/A^ corrections to the matrix element of Qq to be small. 

To perform however a detailed numerical analysis in view of the ratio e' /e, we have to 
include in addition the logarithmic terms which arise from the explicit loop corrections. 
These terms can be calculated in a straightforward way only for the non-factorizable dia- 
grams. To determine the logarithmic behaviour in the factorizable part one has to introduce 
a clear prescription to fix the virtual meson momenta. As we have seen this is not possi- 
ble from a direct matching with the short-distance contributions. This issue is still to be 
investigated. 
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4 Operator Evolution 



To gain a deeper understanding of long-distance effects we may study the corrections 
to density-density operators themselves rather than those to the matrix elements. For that 
we may apply the background field method as used by Fatelo and Gerard in the case 
of current-current operators This approach is powerful as 

• the one-loop effects are given in terms of operators, too. Therefore the analysis is 
process independent (respecting the various decay channels of ^ tttt). 

• the calculation keeps track of the chiral structure in the loop corrections. 

Whereas the method is particularly powerful in the chiral limit on account of the (relatively) 
small number of diagrams which have to be calculated, the inclusion of mass terms is 
difficult, wherefore we restrict our analysis to the case of massless quarks, i.e., to the 
quadratic cutoff dependence. 

Starting from the chiral limit (m^ = 0), in the lagrangian we have to keep a mass term 
for the singlet pseudoscalar 770 arising through the explicit breaking of the Ua{^) symmetry, 

C = :L(Tr(9^[/9^[/t) + !^[Tr(lnf/-lnf/t)]2 
o ^ 12 

= £ + \{d,ed>^n + ^Tr([a^e, ^]d^uu^) - ^miee + oie) . (u) 



To obtain the second line of eq. ([T7|) we decomposed the matrix U in the quantum field ^ 
and the classical field If, 

U = exp{2iC/f) U, e = W , (18) 

U satisfying the equation of motion 

2 

— — — — TD — — — 

Ud'^U^ - d^UU^ = ^Tr(ln f/ - In f/^) ■ 1 , U = exp(2i7r'^A„//) . (19) 

If in the same way we expand the meson density, given in eq. (^ , around the classical field 
(with U ^ W for L ^ R), which yields 

(Dk),, = (Dnh + '-^{W^-j^id'W^ + 2d,U^d'^^ + Wd'0)^^+'^[u^e (20) 



we may calculate the 1/A^ correction to a density-density operator of the form {Dji)ij{DL)ki 
by integrating out the ^ field, the terms linear in ^ in eq. ( pO]) contributing to the non- 
factorizable one-loop diagrams, the quadratic to the factorizable ones, where the rescat- 
tering of the mesons due to the strong interaction may be included using the lagrangian 
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(p!7|). The external legs of a diagram (see fig. 4) now symbolize weak and strong operators, 
respectively, given in terms of the field U. Consequently, the loop corrections preserve the 
chiral structure. 




+ (g) (g) + 






+ <S> (g> — + 





+ — ^ <S) — + 



Fig. 4: One-loop 0{p^) contribution to the long-distance evolution of the operator {Dji)ij{DL)ki, 
the cross circles denoting the densities, the black circles the strong interaction. 

The non-factorizable part (first row of fig. 4) may again be calculated in a straightforward 
(although lengthy) way associating the cutoff to the effective color singlet gauge boson 
connecting the two densities. The factorizable contributions however (second and third 
row of fig. 4), involving trilinear divergences through a certain class of diagrams, show a 
dependence on a shift of the loop momenta {q ^ q + p) already in the quadratic cutoff 
behaviour. This feature, occuring in the study of the operator evolution by means of the 
background field approach, does not appear in the explicit analysis of the various matrix 
elements (see section]^). 

In order to avoid an ambiguity of the operator evolution in the chiral limit, it is sufficient 
to adopt a prescription for the distribution of momenta in the concerned diagrams. Indeed 
we have strong indications for a particular choice (based on arguments beyond the scope 
of this talk), which removes any further shift dependence from the quadratic terms. This 
choice is natural but not necessarily unique. 

Thus performing the analysis, from the sum of the factorizable and the non-factorizable 
parts we obtain the long-distance evolution of the operator Qq, 

4 ^2 _ ^2 



4 Al 



QMlut) = -^^{d,Ud^U%M 1 + 3-fl- + 0{l/N') , (21) 



(4vr/)^ 



the chiral loop corrections being identical to those found in the study of the matrix ele- 
ment (using the background field approach the renormalization factors for the bare decay 
constant / and the bare meson field vr" cancel in the matrix U). 
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Note that the result is given in terms of {dJJd'^tJ^)ds, the latter being the only pseu- 
doscalar operator available at the 0{p^) to describe both, the corresponding current-current 
and density-density four-fermion operators {Q4 and Qq, respectively). 

5 Summary and Conclusion 

We studied the one-loop long-distance contributions to K ^ tttt decay amplitudes apply- 
ing the 1/N expansion in which we used the momentum flow between the currents and 
densities, respectively, to fix the virtual meson momenta. Thus a proper matching with 
the short- distance part is feasible for non-factorizable diagrams. 

The matrix elements of density-density operators, containing factorizable contributions 
which have to be considered explicitly, were found to be clearly determined regarding the 
quadratic cutoff behaviour, yet leaving us with an ambiguity in the logarithmic one. In the 
evolution of the operators themselves, which we investigated by means of the background 
field approach, the corresponding ambiguity already appeared in the quadratic divergences. 
However, using a prescription motivated by theoretical arguments the methods are in 
agreement. 

As a first analytic result we presented the 1/N corrections to the gluon penguin operator 
Qq as well as to the corresponding matrix elements in the chiral limit. The corrections 
are small, provided the change in the value of the hadronic scale is considered. A 
detailed numerical analysis in view of the ratio e'/e requires in addition the inclusion of 
logarithmically divergent terms. This issue is still under investigation. 
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